Decoding triancestral origins, archaic introgression, and natural selection in the Japanese population by whole-genome sequencing.

We generated Japanese Encyclopedia of Whole-Genome/Exome Sequencing Library (JEWEL), a high-depth whole-genome sequencing dataset comprising 3256 individuals from across Japan. Analysis of JEWEL revealed genetic characteristics of the Japanese population that were not discernible using microarray data. First, rare variant-based analysis revealed an unprecedented fine-scale genetic structure. Together with population genetics analysis, the present-day Japanese can be decomposed into three ancestral components. Second, we identified unreported loss-of-function (LoF) variants and observed that for specific genes, LoF variants appeared to be restricted to a more limited set of transcripts than would be expected by chance, with PTPRD as a notable example. Third, we identified 44 archaic segments linked to complex traits, including a Denisovan-derived segment at NKX6-1 associated with type 2 diabetes. Most of these segments are specific to East Asians. Fourth, we identified candidate genetic loci under recent natural selection. Overall, our work provided insights into genetic characteristics of the Japanese population.

Identification of a de novo PUF60 variant associated with craniofacial microsomia.

Craniofacial microsomia (CFM), also known as the oculo-auriculo-vertebral spectrum, is a congenital disorder characterized by hypoplasia of the mandible and external ear due to tissue malformations originating from the first and second branchial arches. However, distinguishing it from other syndromes of branchial arch abnormalities is difficult, and causal variants remain unidentified in many cases. In this report, we performed an exome sequencing analysis of a Brazilian family with CFM. The proband was a 12-month-old boy with clinical findings consistent with the diagnostic criteria for CFM, including unilateral mandibular hypoplasia, microtia, and external auditory canal abnormalities. A heterozygous de novo nonsense variant (c.713C>G, p.S238*) in PUF60 was identified, which was predicted to be pathogenic in silico. PUF60 has been reported as a causal gene in Verheij syndrome, but not in CFM. Although the boy showed craniofacial abnormalities and developmental delay that overlapped with Verheij syndrome, the facial asymmetry with unilateral hypoplasia of the mandible observed in this case did not match the previously reported phenotypes of PUF60 variants. Our findings expand the phenotypic range of PUF60 variants that cover CFM and Verheij syndrome.

Core planar cell polarity genes VANGL1 and VANGL2 in predisposition to congenital vertebral malformations.

Congenital scoliosis (CS), affecting approximately 0.5 to 1 in 1,000 live births, is commonly caused by congenital vertebral malformations (CVMs) arising from aberrant somitogenesis or somite differentiation. While Wnt/ß-catenin signaling has been implicated in somite development, the function of Wnt/planar cell polarity (Wnt/PCP) signaling in this process remains unclear. Here, we investigated the role of Vangl1 and Vangl2 in vertebral development and found that their deletion causes vertebral anomalies resembling human CVMs. Analysis of exome sequencing data from multiethnic CS patients revealed a number of rare and deleterious variants in VANGL1 and VANGL2, many of which exhibited loss-of-function and dominant-negative effects. Zebrafish models confirmed the pathogenicity of these variants. Furthermore, we found that Vangl1 knock-in (p.R258H) mice exhibited vertebral malformations in a Vangl gene dose- and environment-dependent manner. Our findings highlight critical roles for PCP signaling in vertebral development and predisposition to CVMs in CS patients, providing insights into the molecular mechanisms underlying this disorder.

Dyssegmental dysplasia Rolland-Desbuquois type is caused by pathogenic variants in HSPG2 - a founder haplotype shared in five patients.

Dyssegmental dysplasia (DD) is a severe skeletal dysplasia comprised of two subtypes: lethal Silverman-Handmaker type (DDSH) and nonlethal Rolland-Desbuquois type (DDRD). DDSH is caused by biallelic pathogenic variants in HSPG2 encoding perlecan, whereas the genetic cause of DDRD remains undetermined. Schwartz-Jampel syndrome (SJS) is also caused by biallelic pathogenic variants in HSPG2 and is an allelic disorder of DDSH. In SJS and DDSH, 44 and 8 pathogenic variants have been reported in HSPG2, respectively. Here, we report that five patients with DDRD carried four pathogenic variants in HSPG2: c.9970 G > A (p.G3324R), c.559 C > T (p.R187X), c7006 + 1 G > A, and c.11562 + 2 T > G. Two patients were homozygous for p.G3324R, and three patients were heterozygous for p.G3324R. Haplotype analysis revealed a founder haplotype spanning 85,973 bp shared in the five patients. SJS, DDRD, and DDSH are allelic disorders with pathogenic variants in HSPG2.

Association of genetic variation in COL11A1 with adolescent idiopathic scoliosis.

Adolescent idiopathic scoliosis (AIS) is a common and progressive spinal deformity in children that exhibits striking sexual dimorphism, with girls at more than fivefold greater risk of severe disease compared to boys. Despite its medical impact, the molecular mechanisms that drive AIS are largely unknown. We previously defined a female-specific AIS genetic risk locus in an enhancer near the PAX1 gene. Here, we sought to define the roles of PAX1 and newly identified AIS-associated genes in the developmental mechanism of AIS. In a genetic study of 10,519 individuals with AIS and 93,238 unaffected controls, significant association was identified with a variant in COL11A1 encoding collagen (α1) XI (rs3753841; NM_080629.2_c.4004C>T; p.(Pro1335Leu); p=7.07E-11, OR = 1.118). Using CRISPR mutagenesis we generated Pax1 knockout mice (Pax1-/-). In postnatal spines we found that PAX1 and collagen (α1) XI protein both localize within the intervertebral disc-vertebral junction region encompassing the growth plate, with less collagen (α1) XI detected in Pax1-/- spines compared to wild-type. By genetic targeting we found that wild-type Col11a1 expression in costal chondrocytes suppresses expression of Pax1 and of Mmp3, encoding the matrix metalloproteinase 3 enzyme implicated in matrix remodeling. However, the latter suppression was abrogated in the presence of the AIS-associated COL11A1P1335L mutant. Further, we found that either knockdown of the estrogen receptor gene Esr2 or tamoxifen treatment significantly altered Col11a1 and Mmp3 expression in chondrocytes. We propose a new molecular model of AIS pathogenesis wherein genetic variation and estrogen signaling increase disease susceptibility by altering a PAX1-COL11a1-MMP3 signaling axis in spinal chondrocytes.

Adolescent idiopathic scoliosis (AIS) is a twisting deformity of the spine that occurs during periods of rapid growth in children worldwide. Children with severe cases of AIS require surgery to stop it from getting worse, presenting a significant financial burden to health systems and families. Although AIS is known to cluster in families, its genetic causes and its inheritance pattern have remained elusive. Additionally, AIS is known to be more prevalent in females, a bias that has not been explained. Advances in techniques to study the genetics underlying diseases have revealed that certain variations that increase the risk of AIS affect cartilage and connective tissue. In humans, one such variation is near a gene called Pax1, and it is female-specific. The extracellular matrix is a network of proteins and other molecules in the space between cells that help connect tissues together, and it is particularly important in cartilage and other connective tissues. One of the main components of the extracellular matrix is collagen. Yu, Kanshour, Ushiki et al. hypothesized that changes in the extracellular matrix could affect the cartilage and connective tissues of the spine, leading to AIS. To show this, the scientists screened over 100,000 individuals and found that AIS is associated with variants in two genes coding for extracellular matrix proteins. One of these variants was found in a gene called Col11a1, which codes for one of the proteins that makes up collagen. To understand the relationship between Pax1 and Col11a1, Yu, Kanshour, Ushiki et al. genetically modified mice so that they would lack the Pax1 gene. In these mice, the activation of Col11a1 was reduced in the mouse spine. They also found that the form of Col11a1 associated with AIS could not suppress the activation of a gene called Mmp3 in mouse cartilage cells as effectively as unmutated Col11a1. Going one step further, the researchers found that lowering the levels of an estrogen receptor altered the activation patterns of Pax1, Col11a1, and Mmp3 in mouse cartilage cells. These findings suggest a possible mechanism for AIS, particularly in females. The findings of Yu, Kanshour, Ushiki et al. highlight that cartilage cells in the spine are particularly relevant in AIS. The results also point to specific molecules within the extracellular matrix as important for maintaining proper alignment in the spine when children are growing rapidly. This information may guide future therapies aimed at maintaining healthy spinal cells in adolescent children, particularly girls.

Impaired glycine neurotransmission causes adolescent idiopathic scoliosis.

Adolescent idiopathic scoliosis (AIS) is the most common form of spinal deformity, affecting millions of adolescents worldwide, but it lacks a defined theory of etiopathogenesis. Because of this, treatment of AIS is limited to bracing and/or invasive surgery after onset. Preonset diagnosis or preventive treatment remains unavailable. Here, we performed a genetic analysis of a large multicenter AIS cohort and identified disease-causing and predisposing variants of SLC6A9 in multigeneration families, trios, and sporadic patients. Variants of SLC6A9, which encodes glycine transporter 1 (GLYT1), reduced glycine-uptake activity in cells, leading to increased extracellular glycine levels and aberrant glycinergic neurotransmission. Slc6a9 mutant zebrafish exhibited discoordination of spinal neural activities and pronounced lateral spinal curvature, a phenotype resembling human patients. The penetrance and severity of curvature were sensitive to the dosage of functional glyt1. Administration of a glycine receptor antagonist or a clinically used glycine neutralizer (sodium benzoate) partially rescued the phenotype. Our results indicate a neuropathic origin for "idiopathic" scoliosis, involving the dysfunction of synaptic neurotransmission and central pattern generators (CPGs), potentially a common cause of AIS. Our work further suggests avenues for early diagnosis and intervention of AIS in preadolescents.

Family History of Developmental Dysplasia of the Hip is a Risk Factor for the Progression of Hip Osteoarthritis.

BACKGROUND: Developmental dysplasia of the hip (DDH) is considered to have genetic predisposition and presents many intrafamilial occurrences. However, there is no report that evaluates the effect of DDH family history on the progression after the onset of hip osteoarthritis (OA). METHODS: Medical interviews about detailed clinical information including family history were conducted on 298 consecutive patients who had undergone surgery for OA due to DDH. Clinical or radiographic items that are associated with the severity of DDH (total hip arthroplasty [THA], involvement of bilateral DDH, onset age of hip pain, and three radiological indices of DDH: center-edge angle, sharp angle, and acetabular roof obliquity) were collected and evaluated in multivariate analyses for their associations with DDH family history in a qualitative or quantitative manner. Survival time analyses for THA as the endpoint was also performed to evaluate the effects of DDH family history on the progression of OA. RESULTS: The DDH family history showed significant associations with bilateral involvement of DDH (odds ratio = 2.09 [95% confidence interval {CI} 1.05 to 4.16]; P = .037), early onset of hip pain (P = .0065), and radiological severity of DDH (P = .016). The DDH family history showed a significant association with undergoing THA (odds ratio = 2.25 [95% CI 1.09 to 4.66]; P = .029), further supported by the Cox regression analyses (hazards ratio = 1.56 [95% CI 1.15 to 2.11]; P = .0044). CONCLUSION: A DDH family history is a risk factor for the progression of hip OA. Stronger genetic predisposition to DDH leads to faster onset and progression of hip OA.

Compound heterozygous variants in RAB34 in a rare skeletal ciliopathy syndrome.

Skeletal ciliopathies are a heterogenous group of congenital disorders characterized by multiple internal abnormalities, and distinct radiographic presentation. Pathogenic variants in at least 30 cilia genes are known to cause skeletal ciliopathies. Here we report a fetus with an atypical skeletal ciliopathy phenotype and compound heterozygous variants in the RAB34 gene. The affected fetus had multiple malformations, including posterior neck edema, micrognathia, low-set and small ears, auricular hypoplasia, cleft lip and palate, short extremities, and a combination of rarely occurring pre- and postaxial polydactyly. Genome sequencing identified compound heterozygous variants in the RAB34 gene: maternal c.254T>C, p.(Ile85Thr), and paternal c.691C>T, p.(Arg231*) variants. Only the paternal variant was present in the unaffected sibling. Evidence in the literature indicated that Rab34-/- mice displayed a ciliopathy phenotype with cleft palate and polydactyly. These features were consistent with malformations detected in our patient supporting the pathogenicity of the identified RAB34 variants. Overall, this case report further expands genetic landscape of human ciliopathy syndromes and suggests RAB34 as a candidate gene for skeletal ciliopathies.

Interrogating Causal Effects of Body Composition and Puberty-Related Risk Factors on Adolescent Idiopathic Scoliosis: A Two-Sample Mendelian Randomization Study.

Adolescent idiopathic scoliosis (AIS) is the most common form of pediatric musculoskeletal disorder. Observational studies have pointed to several risk factors for AIS, but almost no evidence exists to support their causal association with AIS. Here, we applied Mendelian randomization (MR), known to limit bias from confounding and reverse causation, to investigate causal associations between body composition and puberty-related exposures and AIS risk in Europeans and Asians. For our two-sample MR studies, we used single nucleotide polymorphisms (SNPs) associated with body mass index (BMI), waist-hip ratio, lean mass, childhood obesity, bone mineral density (BMD), 25-hydroxyvitamin D (25OHD), age at menarche, and pubertal growth in large European genome-wide association studies (GWAS), and with adult osteoporosis risk and age of menarche in Biobank Japan. We extracted estimates of the aforementioned SNPs on AIS risk from the European or Asian subsets of the largest multiancestry AIS GWAS (N = 7956 cases/88,459 controls). The results of our inverse variance-weighted (IVW) MR estimates suggest no causal association between the aforementioned risk factors and risk of AIS. Pleiotropy-sensitive MR methods yielded similar results. However, restricting our analysis to European females with AIS, we observed a causal association between estimated BMD and the risk of AIS (IVW odds ratio for AIS = 0.1, 95% confidence interval 0.01 to 0.7, p = 0.02 per SD increase in estimated BMD), but this association was no longer significant after adjusting for BMI, body fat mass, and 25OHD and remained significant after adjusting for age at menarche in multivariable MR. In conclusion, we demonstrated a protective causal effect of BMD on AIS risk in females of European ancestry, but this effect was modified by BMI, body fat mass, and 25OHD levels. Future MR studies using larger AIS GWAS are needed to investigate small effects of the aforementioned exposures on AIS. © 2023 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

Androgen receptor binding sites enabling genetic prediction of mortality due to prostate cancer in cancer-free subjects.

Prostate cancer (PrCa) is the second most common cancer worldwide in males. While strongly warranted, the prediction of mortality risk due to PrCa, especially before its development, is challenging. Here, we address this issue by maximizing the statistical power of genetic data with multi-ancestry meta-analysis and focusing on binding sites of the androgen receptor (AR), which has a critical role in PrCa. Taking advantage of large Japanese samples ever, a multi-ancestry meta-analysis comprising more than 300,000 subjects in total identifies 9 unreported loci including ZFHX3, a tumor suppressor gene, and successfully narrows down the statistically finemapped variants compared to European-only studies, and these variants strongly enrich in AR binding sites. A polygenic risk scores (PRS) analysis restricting to statistically finemapped variants in AR binding sites shows among cancer-free subjects, individuals with a PRS in the top 10% have a strongly higher risk of the future death of PrCa (HR: 5.57, P = 4.2 × 10-10). Our findings demonstrate the potential utility of leveraging large-scale genetic data and advanced analytical methods in predicting the mortality of PrCa.

Evidence of causality of low body mass index on risk of adolescent idiopathic scoliosis: a Mendelian randomization study.

Introduction: Adolescent idiopathic scoliosis (AIS) is a disorder with a three-dimensional spinal deformity and is a common disease affecting 1-5% of adolescents. AIS is also known as a complex disease involved in environmental and genetic factors. A relation between AIS and body mass index (BMI) has been epidemiologically and genetically suggested. However, the causal relationship between AIS and BMI remains to be elucidated. Material and methods: Mendelian randomization (MR) analysis was performed using summary statistics from genome-wide association studies (GWASs) of AIS (Japanese cohort, 5,327 cases, 73,884 controls; US cohort: 1,468 cases, 20,158 controls) and BMI (Biobank Japan: 173430 individual; meta-analysis of genetic investigation of anthropometric traits and UK Biobank: 806334 individuals; European Children cohort: 39620 individuals; Population Architecture using Genomics and Epidemiology: 49335 individuals). In MR analyses evaluating the effect of BMI on AIS, the association between BMI and AIS summary statistics was evaluated using the inverse-variance weighted (IVW) method, weighted median method, and Egger regression (MR-Egger) methods in Japanese. Results: Significant causality of genetically decreased BMI on risk of AIS was estimated: IVW method (Estimate (beta) [SE] = -0.56 [0.16], p = 1.8 × 10-3), weighted median method (beta = -0.56 [0.18], p = 8.5 × 10-3) and MR-Egger method (beta = -1.50 [0.43], p = 4.7 × 10-3), respectively. Consistent results were also observed when using the US AIS summary statistic in three MR methods; however, no significant causality was observed when evaluating the effect of AIS on BMI. Conclusions: Our Mendelian randomization analysis using large studies of AIS and GWAS for BMI summary statistics revealed that genetic variants contributing to low BMI have a causal effect on the onset of AIS. This result was consistent with those of epidemiological studies and would contribute to the early detection of AIS.

Association of genetic variation in COL11A1 with adolescent idiopathic scoliosis.

Adolescent idiopathic scoliosis (AIS) is a common and progressive spinal deformity in children that exhibits striking sexual dimorphism, with girls at more than five-fold greater risk of severe disease compared to boys. Despite its medical impact, the molecular mechanisms that drive AIS are largely unknown. We previously defined a female-specific AIS genetic risk locus in an enhancer near the PAX1 gene. Here we sought to define the roles of PAX1 and newly-identified AIS-associated genes in the developmental mechanism of AIS. In a genetic study of 10,519 individuals with AIS and 93,238 unaffected controls, significant association was identified with a variant in COL11A1 encoding collagen (α1) XI (rs3753841; NM_080629.2_c.4004C>T; p.(Pro1335Leu); P=7.07e-11, OR=1.118). Using CRISPR mutagenesis we generated Pax1 knockout mice (Pax1-/-). In postnatal spines we found that PAX1 and collagen (α1) XI protein both localize within the intervertebral disc (IVD)-vertebral junction region encompassing the growth plate, with less collagen (α1) XI detected in Pax1-/- spines compared to wildtype. By genetic targeting we found that wildtype Col11a1 expression in costal chondrocytes suppresses expression of Pax1 and of Mmp3, encoding the matrix metalloproteinase 3 enzyme implicated in matrix remodeling. However, this suppression was abrogated in the presence of the AIS-associated COL11A1P1335L mutant. Further, we found that either knockdown of the estrogen receptor gene Esr2, or tamoxifen treatment, significantly altered Col11a1 and Mmp3 expression in chondrocytes. We propose a new molecular model of AIS pathogenesis wherein genetic variation and estrogen signaling increase disease susceptibility by altering a Pax1-Col11a1-Mmp3 signaling axis in spinal chondrocytes.

Improved genetic prediction of the risk of knee osteoarthritis using the risk factor-based polygenic score.

BACKGROUND: Polygenic risk score (PRS) analysis is used to predict disease risk. Although PRS has been shown to have great potential in improving clinical care, PRS accuracy assessment has been mainly focused on European ancestry. This study aimed to develop an accurate genetic risk score for knee osteoarthritis (OA) using a multi-population PRS and leveraging a multi-trait PRS in the Japanese population. METHODS: We calculated PRS using PRS-CS-auto, derived from genome-wide association study (GWAS) summary statistics for knee OA in the Japanese population (same ancestry) and multi-population. We further identified risk factor traits for which PRS could predict knee OA and subsequently developed an integrated PRS based on multi-trait analysis of GWAS (MTAG), including genetically correlated risk traits. PRS performance was evaluated in participants of the Nagahama cohort study who underwent radiographic evaluation of the knees (n = 3,279). PRSs were incorporated into knee OA integrated risk models along with clinical risk factors. RESULTS: A total of 2,852 genotyped individuals were included in the PRS analysis. The PRS based on Japanese knee OA GWAS was not associated with knee OA (p = 0.228). In contrast, PRS based on multi-population knee OA GWAS showed a significant association with knee OA (p = 6.7 × 10-5, odds ratio (OR) per standard deviation = 1.19), whereas PRS based on MTAG of multi-population knee OA, along with risk factor traits such as body mass index GWAS, displayed an even stronger association with knee OA (p = 5.4 × 10-7, OR = 1.24). Incorporating this PRS into traditional risk factors improved the predictive ability of knee OA (area under the curve, 74.4% to 74.7%; p = 0.029). CONCLUSIONS: This study showed that multi-trait PRS based on MTAG, combined with traditional risk factors, and using large sample size multi-population GWAS, significantly improved predictive accuracy for knee OA in the Japanese population, even when the sample size of GWAS of the same ancestry was small. To the best of our knowledge, this is the first study to show a statistically significant association between the PRS and knee OA in a non-European population. TRIAL REGISTRATION: No. C278.

Microtubule stabilization targeting regenerative chondrocyte cluster for cartilage regeneration.

Purpose: Chondrocytes (CHs) in cartilage undergo several detrimental events during the development of osteoarthritis (OA). However, the mechanism underlying CHs regeneration involved in pathogenesis is largely unknown. The aim of this study was to explore the underlying mechanism of regeneration of CHs involved in the pathological condition and the potential therapeutic strategies of cartilage repair. Methods and Materials: CHs were isolated from human cartilage in different OA stages and the high-resolution cellular architecture of human osteoarthritis was examined by applying single-cell RNA sequencing. The analysis of gene differential expression and gene set enrichment was utilized to reveal the relationship of cartilage regeneration and microtubule stabilization. Microtubule destabilizer (nocodazole) and microtubule stabilizer (docetaxel) treated-human primary CHs and rats cartilage defect model were used to investing the effects and downstream signaling pathway of microtubule stabilization on cartilage regeneration. Results: CHs subpopulations were identified on the basis of their gene markers and the data indicated an imbalance caused by an increase in the degeneration and disruption of CHs regeneration in OA samples. Interestingly, the CHs subpopulation namely CHI3L1+ CHs, was characterized by the cell regenerative capacity, stem cell potency and the activated microtubule (MT) process. Furthermore, the data indicated that MT stabilization was effective in promoting cartilage regeneration in rats with cartilage injury model by inhibiting YAP activity. Conclusion: These findings lead to a new understanding of CHs regeneration in the OA pathophysiology context and suggest that MT stabilization is a promising therapeutic target for OA and cartilage injury.

Null and missense mutations of ERI1 cause a recessive phenotypic dichotomy in humans.

ERI1 is a 3'-to-5' exoribonuclease involved in RNA metabolic pathways including 5.8S rRNA processing and turnover of histone mRNAs. Its biological and medical significance remain unclear. Here, we uncover a phenotypic dichotomy associated with bi-allelic ERI1 variants by reporting eight affected individuals from seven unrelated families. A severe spondyloepimetaphyseal dysplasia (SEMD) was identified in five affected individuals with missense variants but not in those with bi-allelic null variants, who showed mild intellectual disability and digital anomalies. The ERI1 missense variants cause a loss of the exoribonuclease activity, leading to defective trimming of the 5.8S rRNA 3' end and a decreased degradation of replication-dependent histone mRNAs. Affected-individual-derived induced pluripotent stem cells (iPSCs) showed impaired in vitro chondrogenesis with downregulation of genes regulating skeletal patterning. Our study establishes an entity previously unreported in OMIM and provides a model showing a more severe effect of missense alleles than null alleles within recessive genotypes, suggesting a key role of ERI1-mediated RNA metabolism in human skeletal patterning and chondrogenesis.

Genetic factors affecting survival in Japanese patients with sporadic amyotrophic lateral sclerosis: a genome-wide association study and verification in iPSC-derived motor neurons from patients.

BACKGROUND: Several genetic factors are associated with the pathogenesis of sporadic amyotrophic lateral sclerosis (ALS) and its phenotypes, such as disease progression. Here, in this study, we aimed to identify the genes that affect the survival of patients with sporadic ALS. METHODS: We enrolled 1076 Japanese patients with sporadic ALS with imputed genotype data of 7 908 526 variants. We used Cox proportional hazards regression analysis with an additive model adjusted for sex, age at onset and the first two principal components calculated from genotyped data to conduct a genome-wide association study. We further analysed messenger RNA (mRNA) and phenotype expression in motor neurons derived from induced pluripotent stem cells (iPSC-MNs) of patients with ALS. RESULTS: Three novel loci were significantly associated with the survival of patients with sporadic ALS-FGF1 at 5q31.3 (rs11738209, HR=2.36 (95% CI, 1.77 to 3.15), p=4.85×10-9), THSD7A at 7p21.3 (rs2354952, 1.38 (95% CI, 1.24 to 1.55), p=1.61×10-8) and LRP1 at 12q13.3 (rs60565245, 2.18 (95% CI, 1.66 to 2.86), p=2.35×10-8). FGF1 and THSD7A variants were associated with decreased mRNA expression of each gene in iPSC-MNs and reduced in vitro survival of iPSC-MNs obtained from patients with ALS. The iPSC-MN in vitro survival was reduced when the expression of FGF1 and THSD7A was partially disrupted. The rs60565245 was not associated with LRP1 mRNA expression. CONCLUSIONS: We identified three loci associated with the survival of patients with sporadic ALS, decreased mRNA expression of FGF1 and THSD7A and the viability of iPSC-MNs from patients. The iPSC-MN model reflects the association between patient prognosis and genotype and can contribute to target screening and validation for therapeutic intervention.

A novel CCDC91 isoform associated with ossification of the posterior longitudinal ligament of the spine works as a non-coding RNA to regulate osteogenic genes.

Ossification of the posterior longitudinal ligament of the spine (OPLL) is a common intractable disease that causes spinal stenosis and myelopathy. We have previously conducted genome-wide association studies for OPLL and identified 14 significant loci, but their biological implications remain mostly unclear. Here, we examined the 12p11.22 locus and identified a variant in the 5' UTR of a novel isoform of CCDC91 that was associated with OPLL. Using machine learning prediction models, we determined that higher expression of the novel CCDC91 isoform was associated with the G allele of rs35098487. The risk allele of rs35098487 showed higher affinity in the binding of nuclear proteins and transcription activity. Knockdown and overexpression of the CCDC91 isoform in mesenchymal stem cells and MG-63 cells showed paralleled expression of osteogenic genes, including RUNX2, the master transcription factor of osteogenic differentiation. The CCDC91 isoform directly interacted with MIR890, which bound to RUNX2 and decreased RUNX2 expression. Our findings suggest that the CCDC91 isoform acts as a competitive endogenous RNA by sponging MIR890 to increase RUNX2 expression.

Exome Sequencing in Monogenic Forms of Rickets.

OBJECTIVE: To understand the phenotypic and genotypic spectrum of genetic forms of rickets in 10 families. METHODS: Detailed clinical, radiographic, and biochemical evaluation of 10 families with phenotypes suggestive of a genetic cause of rickets was performed. Molecular testing using exome sequencing aided in the diagnosis of six different forms of known genetic causes. RESULTS: Eleven disease-causing variants including five previously reported variants (CYP27B1:c.1319_1325dup, p.(Phe443Profs*24), VDR:c.1171C>T, p.(Arg391Cys), PHEX: c.1586_1586+1del, PHEX: c.1482+5G>C, PHEX: c.58C>T, p.(Arg20*)) and six novel variants (CYP27B1:c.974C>T, p.(Thr325Met), CYP27B1: c.1376G>A, p.(Arg459His), CYP2R1: c.595C>T, p.(Arg199*), CYP2R1:c.1330G>C, p.(Gly444Arg),SLC34A3:c.1336-11_1336-1del, SLC2A2: c.589G>C, p.(Val197Leu)) in the genes known to cause monogenic rickets were identified. CONCLUSION: The authors hereby report a case series of individuals from India with a molecular diagnosis of rickets and provide the literature review which would help in enhancing the clinical and molecular profile for rapid and differential diagnosis of rickets.

Identification of a Functional Susceptibility Variant for Adolescent Idiopathic Scoliosis that Upregulates Early Growth Response 1 (EGR1)-Mediated UNCX Expression.

Adolescent idiopathic scoliosis (AIS) is a serious health problem affecting 3% of live births all over the world. Many loci associated with AIS have been identified by previous genome wide association studies, but their biological implication remains mostly unclear. In this study, we evaluated the AIS-associated variants in the 7p22.3 locus by combining in silico, in vitro, and in vivo analyses. rs78148157 was located in an enhancer of UNCX, a homeobox gene and its risk allele upregulated the UNCX expression. A transcription factor, early growth response 1 (EGR1), transactivated the rs78148157-located enhancer and showed a higher binding affinity for the risk allele of rs78148157. Furthermore, zebrafish larvae with UNCX messenger RNA (mRNA) injection developed body curvature and defective neurogenesis in a dose-dependent manner. rs78148157 confers the genetic susceptibility to AIS by enhancing the EGR1-regulated UNCX expression. © 2022 American Society for Bone and Mineral Research (ASBMR).

SLC4A2, another gene involved in acid-base balancing machinery of osteoclasts, causes osteopetrosis.

SLC4A2 belongs to the Na+-independent solute carrier family 4 (SLC4) of anion exchangers, which regulate electroneutral exchange of Cl- for HCO3- and mediate intra- and extra-cellular pH, chloride concentration and cell volume. Slc4a2 also participates in gastric acid secretion, spermatogenesis and osteoclastogenesis. During osteoclast differentiation, Slc4a2 is exclusively expressed at the contra-lacunar membrane and is up-regulated with osteoclast maturation. Bi-allelic Slc4a2 loss-of-function mutations have been known to cause osteopetrosis in mice and cattle, but not in human. Recently, we have identified bi-allelic pathogenic variants in SLC4A2 in a patient affected by osteopetrosis with severe renal insufficiency, suggesting SLC4A2 deficiency causes a new type of autosomal recessive osteopetrosis (osteopetrosis, Ikegawa type). In this article, we review the advances in exploring the multiple functions of SLC4A2 with emphasis on its roles in osteoclast. Our review would contribute to understanding of the phenotypic spectrum and the pathomechanism of SLC4A2-associated osteopetrosis.